Accuracy of topological entanglement entropy on finite cylinders 
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Topological phases are unique states of matter which support non-local excitations which behave 
as particles with fractional statistics. A universal characterization of gapped topological phases is 
provided by the topological entanglement entropy (TEE). We study the finite size corrections to the 
TEE by focusing on systems with Z2 topological ordered state using density-matrix renormalization 
group (DMRG) and perturbative series expansions. We find that extrapolations of the TEE based 
on the Renyi entropies with Renyi index n > 2 suffer from much larger finite size corrections than 
do extrapolations based on the von Neumann entropy. In particular, when the circumference of the 
cylinder is about ten times the correlation length, the TEE obtained using von Neumann entropy has 
an error of order 10 -3 , while for Renyi entropies it can even exceed 40%. We discuss the relevance 
of these findings to previous and future searches for topological ordered phases, including quantum 
spin liquids. 



Topological phases are exotic states of matter that 
are characterized by ground state degeneracy dependent 
upon global topology of the system on which the phase 
resides, and which host exotic excitations with fractional 
quantum statistics. In recent years, surprising connec- 
tions have emerged between topological phases and quan- 
tum information, stimulated by the prospect of using 
them to construct an inherently fault-tolerant quantum 
computer [lj, |2j. Two dimensional phases with topologi- 
cal order are well known in connection with the fractional 
quantum Hall effect [3(, but are also expected to exist in 
frustrated quantum magnets [J, [5j. 

While topological phases are not characterized by any 
local order parameter, theory shows that they can be 
identified by non-local quantum entanglement, specifi- 
cally the topological entanglement entropy (TEE) of the 
ground states [y, [Z| . The entanglement entropy of a subre- 
gion A of the system with a smooth boundary of length L 
is defined from the reduced density matrix p A , according 
to 

Si(A) = -Tr\p A hi(j> A )], 

and takes the form 

S 1 {A)=a 1 L- 1 , 

where 7 is the TEE. Severe finite-size corrections of 
the formulations in Ref. [g, lZ] due to lattice-scale effects 
greatly hinder their practical application [8|. Instead, two 
of us have recently proposed a practical and extremely 
simple scheme, the "cylinder construction" , to accurately 
calculate TEE[9(. The cylinder construction simply con- 
sists of using the DMRG [10] to calculate the usual von 
Neumann entanglement entropy for the division of a 
cylinder into two equal halves by a flat cut, and ex- 
tracting the TEE from its asymptotic large-circumference 
limit. Thereby, we can practically identify topological 
phases in arbitrary realistic models, including physical 
spin models |9l 1 11|. 



The work above utilized the von Neumann entangle- 
ment entropy, and achieved an accurate extrapolation of 
the TEE term. In the literature, many works study in- 
stead the generalized Renyi entropies, defined as 



S n (A) = 



1 



1 -n 



lnTr^, 



while the von Neumann entropy Si is defined as the limit 
n — > 1. For simplicity, we will call S n as Renyi entangle- 
ment entropy when n > 2, while von Neumann entangle- 
ment entropy when n = 1. Theoretically, the universal 
TEE is expected to obtain also for the Renyi entropy, 
with 7 independent of the Renyi index. However, ex- 
trapolations in the literature based on the Renyi entropy 
appear to be substantially less accurate than those based 
on the von Neumann entropy, even for larger boundary 
lengths L[12|. In particular, the extrapolated TEE in 
Ref.[l2j from the second-order Renyi entropy S2 devi- 
ates from the expected value with an error an order of 
magnitude larger than that from von Neumann entangle- 
ment entropy Si Q . This suggests that extrapolations of 
Renyi entropies have significantly larger finite-size effects 
than von Neumann entropy. 

In this letter, we study the finite-size effects in the TEE 
systematically for two canonical models of phases with Z2 
topological order, and confirm the above suggestion. We 
attempt to cast our results in terms of the expected form, 



S n (L) = a n L - 7„ 



(1) 



as a function of Renyi index n. First, we study the Toric- 
Code model whose TEE is known, using DMRG and per- 
turbative series expansions 13] . We then turn to the more 
realistic S — 1/2 antiferromagnetic Heisenberg model on 
the Kagome lattice. For both cases, we find that the 
Renyi entropies do have substantially larger finite-size 
corrections than the von Neumann entropy. We provide 
some understanding of this tendency from the fact, which 
we show from the series expansion, that the line term a n 
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FIG. 1: (Color online) The extrapolated TEE 7™ using von 
Neumann entropy (i.e., n = 1) and Renyi entropy (i.e., n > 
2) for the Toric-Code model, as a function of the applied 
magnetic field h x , for (a) h z — and (b) symmetric case 
h z — h x . The dashed lines represent the expected universal 
value, In 2, for the TEE in the thermodynamic limit. 



FIG. 2: (Color online) Moving two point data fits for the ex- 
trapolated TEE 7„ for the Toric-Code model for (a) h x — 0.3 
with hz = 0, and (b) h x = h z — 0.3, as a function of the cylin- 
der circumference L y and entropy index n. Here, 'fn(Ly) are 
fitted using two data points, i.e., S„(L y ) and S n (L y + 2) using 
Eq.JTJ. The dashed lines represent the expected asymptotic 
value 7 == In 2. 



varies more rapidly with parameters with increasing n. 
This makes the extraction of the sub-dominant -f n term 
less reliable. 

Toric-Code Model We begin with the well known 
Toric-Code model (TCM)[1], (See Supplementary Infor- 
mation for notational details). Without applied fields, 
i.e., h x = h z = 0, the pure TCM is exactly solublefl] with 
a Z 2 topological ordered ground state. After turning on 
the magnetic fields, the model is no longer exactly solu- 
ble. Previous studies 14H16| show that the Z 2 topolog- 



ical phase remains stable and robust until the magnetic 
fields are large enough that the system undergoes a phase 
transition from the topological phase to a trivial one. Nu- 
merically, Jiang, Wang and Balents have systematically 
calculated the von Neumann entanglement entropy Si 
using cylinder construction [9j, and extrapolated an ac- 
curate TEE 71 = ln(2) in the Z 2 topological phase, even 
very close to the phase transition point. In this paper, we 
further calculate the Renyi entropy to study the finite- 
size effects on the TEE j n as a function of Renyi index 
n. To make sure that S n only scales with the cylinder 
circumference L y — L, we will work in the long cylinder 
limit, i.e., cylinder length L x — > 00. Therefore, we can 
directly extrapolate the TEE j n using Eq.([T|). 

To see the finite-size corrections, we first consider an 
applied magnetic field only along x-direction, i.e., h z — 0. 
The extrapolated TEE j n are shown in FigQ] (a) as a 
function of magnetic field h x for fields within the topo- 
logical phase. As shown in the Supplementary Material, 
for the fields in Fig. [TJ the spin correlation length remains 
of order one lattice spacing or smaller. We see that the 
linear fit using data for L y = 6 ~ 12 and Eq.([T]) give s 
quite accurate results for 71 for h x < h x rs 0.328[14j. 
Even for h x — 0.3, very close to the quantum phase tran- 
sition, we obtain 71 = 0.6945(20), which is accurate to 
the expected value (i.e., the dashed line) to a fraction 
of percent ~ 0.2%. By contrast, the estimated TEE j n 



obtained from the Renyi entropy for n > 2 shows dramat- 
ically larger deviations from the universal value. These 
deviations grow when approaching the phase transition 
point, where the errors become an order of magnitude 
larger than those for 71. For example, 72 = 0.714(5) at 
h x = 0.3, with an error around ~ 3%. The deviations 
also increase with Renyi index n, e.g., ~ 7% for n = 4, 
as shown in Fig[T] (a) and inset of FigJU Similar results 
hold for the symmetric case h x — h z , with even larger 
finite-size corrections, as shown in FigJIJb) and inset of 
FiglU Systematically, the finite-size corrections to the 
TEE 7„ defined by Eq.{T]) are much larger for the Renyi 
entropy than for the von Neumann entropy. 

Our expectation is that in the thermodynamic limit 
L y — > 00, all 7„ should converge properly to the univer- 
sal value. We look for signs of this tendency, by using a 
moving two data-point fit. We can define "f n (L y ) using 
two data points with different cylinder circumferences L y 
and L y + 2. Examples of such ~f n (L y ) are shown in Fig[2] 
for h x = 0.3 and h z = in (a), and for h z = h x = 0.3 in 
(b), as a function of L y . For both cases, 71 (L y ) quickly 
converges to the expected value, which is in sharp con- 
trast to 7„(L y ) using Renyi entropies with n > 2. For the 
latter cases, the TEE is systematically over-estimated. 
The curves in Fig. [5] do at least show downward curva- 
ture, consistent with eventual convergence to the univer- 
sal value for larger L y , but it is clear the much larger 
systems would be required to test this in detail. 

To understand the origin of the n-dependence in the 
extrapolations, we turn to a perturbative series expan- 
sion calculation of the line term a n . Since the TEE is 
obtained after subtraction of the much larger line term, 
an accurate extrapolation of 7„ also requires an accurate 
calculation of a n . While there have been perturbative 
studies of the Kitaev model in a field 17| and some entan- 
glement properties have also been studied[18|, [l9(, we are 
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weight of a cluster W n (c) is defined recursively by the 
relations 




FIG. 3: (Color online) Line entropy a n obtained from DMRG 
(block square) and perturbative calculation (red circle) for 
both von Neumann entropy (i.e., n = I) and Renyi entropy 
(i.e., n > 2) for the Toric-Code model, as a function of the 
magnetic field h x , here h z = 0. 



not aware of any exact perturbative evaluation of its line 
entropy relevant to our geometry. In order to carry out 
this calculation, we turn to Linked Cluster Methods. 13 1 

We consider one of the ground state sectors of the TCM 
and map the problem on to the Transverse-Field Ising 
Model (TFIM) if only h x is non-zero 14J and on to a Z 2 
Lattice-gauge model if both h x and h z are non-zero. [16j 
In either case, one has a unique non-degenerate ground 
state. Here, we will consider only the former mapping 
and the case with only h x non-zero. Up to 4-th order, the 
dependence on h x and h z are additive and hence know- 
ing the dependence on h x and the symmetry under the 
interchange of h x and h z one can easily write down the 
full dependence on h x and h z . 

Given the exact mapping between the models, if we 
were to calculate some property like the ground state en- 
ergy in a series expansion in the field, the expansions 
would be identical to the TFIM. However, there is a 
crucial difference for the entanglement entropies. |20| The 
TFIM variables sit at the center of the plaquettes whereas 
the TCM variables and the perturbing fields live on the 
bonds. Furthermore, each state in the TCM is a linear 
superposition of many basis states (say in the a z basis). 
This affects how the states are represented on either side 
of the partition and hence the reduced density matrix as 
we discuss below. 

In the linked cluster method, we can define a clus- 
ter by a set of bonds, where the perturbative fields are 
present. 2l|, |22| The entanglement entropies can be ex- 
pressed as 



S n =£y„(c), (2) 

c 

where the sum is over all possible clusters c and the 



W n (c)=S n (c)-J2Wn(s), 



(3) 



where S n (c) is the entanglement entropy for the cluster 
c and the sum over s is over all subclusters of c. 

One can show that, for the line entropy, the only clus- 
ters that will give non-zero contributions in powers of h 
are those that are (i) linked and (ii) that contain at least 
one bond in subsystem A and one in subsystem B.[2lLl22j| 
Here, one should note that two bonds are linked if they 
meet at a site or if they are on the opposite sides of an 
elementary plaquette (because they can both change the 
flux through a common plaquette). This implies that 
all such clusters must be situated close to the interface 
between A and B. Since such linked clusters can be trans- 
lated along the line, it follows that, for a large system, 
this entropy is proportional to the length of the line. In 
fourth order, we can group the perturbations into just 
two distinct clusters, whose calculational details can be 
found in the Supplementary materials. 

For line term of the Renyi entropy for n > 1, to order 
h A , we obtain 



-(In 2 h A 

2 V 32 



3n h A . 
i-l 128 ; 



while, the von-Neumann entropy (n — » 1) becomes 

e*i = i(m2 - (— - — ln2)h 4 - — h 4 In ft). 
2 V l 128 32 ; 32 ' 



(4) 



(5) 



Note that the innocuous h n In h singularity is inevitable 
for the von-Neumann entropy in any model, where there 
are Schmidt states, whose weight is zero in the unper- 
turbed model but becomes non-zero as a power of the 
perturbation parameter. 

In. Fig(3] we show the line entropies a n obtained in 
DMRG compared with the series expansion results (See 
also Fig |S3l in Supplementary materials) . The agreement 
is excellent. The important thing to note is the n de- 
pendence of the line entropy. The linear n dependence 
in Eq.Q means that with increasing n, the line entropy 
changes more rapidly with the applied fields. Since, the 
topological entanglement entropy is obtained after sub- 
traction of the much larger line entropy, it follows that 
with increasing n the entanglement entropy would have 
a much larger finite size correction, as the correlation 
length increases. 

Kagome Heisenberg Model We now turn to the 
spin-1/2 Heisenberg model on the Kagome lattice, for 
which compelling [23J, |2J] and direct evidence [9|, [l2| for 
a Zi topological quantum spin liquid has been obtained 
by extensive DMRG studies. In particular, highly ac- 
curate TEE 71 = ln(2) has been obtained using cylin- 
der construction [9|, for the model with both first- and 




FIG. 4: (Color online) The extrapolated TEE j n for the 
Kagome J1-J2 Heisenberg model as a function of entropy in- 
dex n at both J2 = 0.10 and J2 = 0.15. Inset: 7„ for TCM as 
a function of entropy index n, for the case h x — 0.3, h z = 
(black square), and the symmetric case h x = h z = 0.3 (red 
circle). The dashed lines correspond to ln(2). 



second-neighbor interactions - see the Supplementary 
Information for details of the definition. Specifically, 
the extrapolated TEE 71 = 0.698(8) at J 2 = 0.10 and 
71 = 0.694(6) at J 2 = 0.15, both within one percent of 
ln(2) = 0.693. 

As in the TCM, we find that the Renyi entropies give 
much less accurate estimates of j n . The extrapolated 
TEE 7„ for n > 2 clearly deviates from the expected 
value, even when the cylinder circumference is much 
larger than the correlation length, i.e., L y m 10£ (the 
correlation lengths are known from the earlier study in 
Ref.[9j). For example, as shown in Fig |S51 a linear fit us- 
ing Eq.© gives 72 w 0.42(1) at both J 2 = 0.10 and 
J 2 = 0.15: a huge error of ~ 40%. Moreover, with 
increasing Renyi index n, the deviation becomes even 
larger, reaching, for example ~ 60% for n = 4, as shown 
in Fig|4] These results show that large finite-size correc- 
tions to the Renyi entropies obtain not only in the "ar- 
tificial" TCM, but also in realistic quantum spin Hamil- 
tonians. 

Summary and Conclusion In this paper, we stud- 
ied the finite-size scaling of the TEE for systems with 
Z 2 topological order, using DMRG simulations and per- 
turbative series expansions. We find that generally the 
finite-size errors in the TEE based on the Renyi entropy 
estimators (i.e., n > 2) are much larger than those ob- 
tained from the von Neumann entropy (i.e., n — I). In 
particular, when the cylinder circumference is around ten 
times the correlation length, L y « 10£, the extrapolated 
TEE using von Neumann entropy is quite accurate with 
an error of order 10 -3 . On the contrary, the error can 
be orders of magnitude larger for Renyi entropy. For in- 



stance, for spin-1/2 Kagome Heisenberg model, the error 
is only around a fraction of percent for von Neumann en- 
tropy, while it is 40% or even larger for Renyi entropy. 
Perturbative study of the TCM shows that the larger fi- 
nite size corrections for the TEE originates in part from 
the enhanced variation with parameters of the line en- 
tropy a n with increasing n. This indicates moreover that 
estimates of the TEE become less accurate with increas- 
ing n. We note that errors of the magnitude found here 
for n > 2 in the Kagome Heisenberg model are large 
enough to perhaps preclude a definitive identification of 
the topological phase, even if we assume that a universal 
value obtains in the thermodynamic limit. Our results 
clearly indicate that great care must be taken into ac- 
count for finite size corrections in numerical calculations 
of the TEE, particularly those based on Renyi entropies 
with 71 > 2. This gives techniques, such as DMRG, which 
have direct access to the full density matrix and hence 
von Neumann entropy, a distinct advantage. Although 
in this paper we have focused on systems with Z 2 topo- 
logical order, similar conclusions may be expected more 
generally. 
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I. Toric-Code Model in Magnetic Field 




FIG. SI: The entanglement entropy for the Toric-Code model 



ai L x — oo. Here n is the entropy index. Inset: Square lattice 
with L x — 10 and L y — 6. Here S represents the star operator 
A s , while P represents the plaquette operator B p . 



The Toric-Code modeljlj with an applied magnetic 
field is given by 

H = -J s Y,A s -J p J2B P -h x J2^-h z J2^( S1 ) 



where of and of are Pauli matrices, A s — n ies crf and 
B p = IL; e paf. Subscripts s and p refer respectively to 
vertices and plaquettes on the square lattice, whereas i 
runs over all bonds where spin degrees of freedom are 
located. Without magnetic field, i.e., h x — h z = 0, the 
pure TCM can be solved exactly [l|. It exhibits a 4- fold 
ground state degeneracy on the torus, and has Zi topo- 
logical order with total quantum dimension D = 2. All 
elementary excitations are gapped and characterized by 
eigenvalues A s — — 1 (a Z 2 charge on site s) and B p = — 1 
(a Z2 vortex on plaquette p) . After turning on the mag- 
netic fields, the model cannot be solved exactly anymore. 
However, previous studies [14rll6| show that the Z 2 topo- 
logical phase remains quite stable and robust until the 
magnetic fields are large enough to induce a phase transi- 
tion from the topological phase to the trivial one. Specif- 
ically, such a phase transition takes place at the critical 
magnetic field h x ~ 0.32 when h z = 0, while h x « 0.34 
along the symmetric line with h z = h x . 

For the DMRG simulation, we consider an equivalent 
square lattice, where the spin operators a x and a z sit on 
the sites instead of the bonds. Therefore, the star op- 
erator A s and the plaquette operator B p of the original 
lattice now sit on alternating plaquettes in the equivalent 
square lattice, as shown in Figure IS1[ labeled as S and 
P, respectively. Note that on this equivalent square lat- 
tice, there are an even number of dangling spins within 
each plaquette at the open edges. For the pure Toric- 
Code model with cylinder boundary condition, the first 




FIG. S2: The spin-spin correlation length for the Toric- 
Code model as a function of the applied magnetic field h x in 
Eq. (|Sl[) . with h z = 0.0. Inset: Examples of spin-spin correla- 
tion function ((Toof) = ((Toof) — (< 7 o)(°f) along the cylinder, 
at different magnetic field h x . 



2 L v/ 2 ~ 1 eigenvalues of the reduced density matrix pa are 
degenerate and equal to \/2 Ly / 2 ~ l , while all the other 
eigenvalues are zero. Therefore, the entanglement en- 
tropy, defined as S n = -^j In Trp^ , does not depend on 

the entropy index n, and is equal to S n = ":, ' L y — ln(2), 
with coefficient a n = ln(2)/2 and TEE 7„ = ln(2). 

After turning on the magnetic field, the degeneracy 
of entanglement spectrum is lifted, and the correlation 
length £ becomes finite, as shown in Fig|S2l (note that 
other types of the correlation length may be defined, but 
are expected to differ just by a factor of order 1 from the 
above one) . Therefore, the entanglement entropy S n 
will depend on the entropy index n, and it follows that 



0.01 





FIG. S4: Graphs needed for the perturbative calculation of 
the line entropy to 4th order. The dashed line separates sub- 
systems A and B. The solid lines denote the perturbing fields 
in the cluster. 



der series expansion. The series expansion results match 
the DMRG values remarkably accurately, though as ex- 
pected, the difference between the two increases markedly 
on approaching the phase transition point, especially for 
larger n. 



FIG. S3: Deviation of the coefficient Sa„ = a n (h x = 0) — 
a n (h x ) obtained from DMRG using Eq.fTJ (red circle) and 
theoretical formula in Eq.Q up to 4th order series expansion 
for Renyi entropy (black square) with (a) n — 2, (b) n = 3 
and (c) n = 4 for the Toric-Code model, as a function of the 
applied magnetic field h x . Here h z = 0. 



so too will the coefficient a n and the extrapolated TEE 
7„. Indeed, our results show that when approaching the 
phase transition point (e.g., h x — 0.34 along the symmet- 
ric line with h x = h z ), the dependence of the TEE 7„ on 
n becomes clearer and clearer, as shown in FigJIJb) in 
the main text, and Fig lSll here. For example, the corre- 
lation length £ sa 1 lattice spacing at h x = 0.30, and the 
resulting TEE 71 = 0.696(4) obtained from the von Neu- 
mann entanglement entropy is still quite accurate with 
an error around a few fraction of 10~ 3 , i.e., ~ 0.4%. On 
the contrary, the fitted TEE 7„ using Renyi entropy (i.e., 
n > 2) deviates from the expected value clearly, with an 
error more than order of magnitude larger than for 71. 
For instance, 72 = 0.726(6), with an error around 5%, 
around one order of magnitude larger than that for 71. 
Such a deviation is even larger with the increase of n, e.g., 
74 = 0.762(8), with an error around 10%, more than one 
order of magnitude larger than that for 71 . 

Deviations are also found in the line entropies a n . 
Fig |S3l shows the comparison for a„ obtained numerically 
using Eq.JT]) and theoretically using Eq.Q up to 4th or- 



II. Perturbative calculations 

The perturbative calculations to order h 4 for the Toric 
Code Model require just two clusters shown in Fig|S4l 
(i) A cluster consisting of all four bonds incident on a 
boundary site between A and B and (ii) A cluster with 
two bonds across a plaquette facing each other across 
the boundary between A and B. The first has a count 
of one per boundary site. The second has a count of 2 
per boundary site. Once grouped this way, no subgraph 
subtraction is needed in 4th order. 

We now discuss the calculation of S n (c), the Renyi en- 
tropy when only the fields on the bonds of a cluster c are 
present. [22j First the ground state wavefunction is calcu- 
lated perturbatively in the TFIM variables, with all the 
perturbing fields present in the cluster. Then, for every 
basis state in the TFIM model, a representative a z basis 
state is obtained in the TCM variable. When all the 4 
perturbative fields are present at a site, the states in the 
TCM must be enlarged by the operation of the projec- 
tion operator 1 + A s at the site. This will give us the 
ground state for the cluster in the TCM variables. Once, 
the ground state for the cluster is known in the TCM 
variables, it is straighforward to calculate the reduced 
density matrix for subsystem A and hence the entangle- 
ment entropies, S n (c). 

First we consider graph (ii) from the figure. This graph 
has just two bonds, one in A and one in B. In the Ian- 



guage of TFIM, there are three plaqucttes on which Vi- 
sons can be created or destroyed. However, since Vi- 
sons are created or destroyed in pairs, there are only 4 
relevant states in the Hilbert space. These can be de- 
noted by the number of Visons present in the plaquette 
as |1 >= |0, 0,0 >, |2 >= |1, 1,0 >, |3 >= |0, 1, 1 > 
and 1 4 >= |1, 0, 1 >. These four states can be related to 
the TCM in terms of o z on the two bonds of the cluster. 
The four states become |1 >= |1, 1 >, |2 >= | — 1, 1 >, 
|3 >= 1 1, — 1 > and |4 >= | - 1,-1 >. In this 4- 
dimcnsional Hilbert space, the Hamiltonian matrix is 



H = 



-h -h 

-h 4 -h 

-h 4 -h 

-h -h 4 



The ground state wavefunction correct to order required 
for h A evaluation of Renyi entropies is 
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Let Ai label the basis states for subsystem A and Bk 
label basis states for subsystem B. Then reduced density 
matrix for subsystem A has matrix elements 



(Pa) 



Ai.Ai 



Bk 



i>A i ,B k i>A i 



(S2) 



In the TCM representation, there are just two states 
on the A side, and the reduced density matrix becomes 



PA 



1 h 2 h 4 

1 16 128 
h 

4 



h 

4 
h 2 , h 4 

16 ' 128 



From this it follows that the two eigenvalues of the 

,4 > 4 

reduced density matrix are 1 — ^g and ^gg. It leads to 
Renyi entanglement entropy of 
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—nh A /i 4 
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Now, we turn to graph (i) in figure, which has 4 bonds 
all sharing a boundary site between A and B. There are 4 
plaqucttes on which Visons can be created or destroyed. 
But, because they are created or destroyed in pairs, there 
number is conserved modulo 2. Thus, there are 8 relevant 
basis states. We can order the plaquettes in a clock- 
wise manner. Then the 8 states are: |1 >= |0, 0,0,0 >, 
|2>=|1,1,0,0>, |3>= |0,1,1,0 >, |4>=|0,0,1,1>, 
|5 >= |1, 0,0,1 >, |6 >= |1, 1,1,1 >, |7 >= |1, 0,1,0 > 
and |8 >= |0, 1,0, 1 >. The 8-dimensional Hamiltonian 



matrix is given by 
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It leads to ground state wave function needed for 4th 
order calculation of Renyi entropies of: 
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Transforming to the TCM variables a z , there are 4 states 
on the A side, and remembering that each TFIM state is 
a superposition of two TCM basis states after applying 
the projection operator (1+A S ), the 4 x 4 reduced density 
matrix becomes: 
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From this it follows that the eigenvalues of 2^ are 1 



3hT__9h^ 1 j. 3h_ J 7fe 

4 64 ' ~r 4 ~r 64 ' 

entanglement entropies: 
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and |j . These lead to Renyi 
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Note that this graph also gives the full line entropy of 
the unperturbed TCM, which is In 2 for every site that 
is shared between subsystems A and B. For the unper- 
turbed model, it is easily shown that any larger graph 
will give zero upon subgraph subtraction as no additional 
degeneracy arises on the line. 



III. Kagome Heisenberg Model 

The spin-1/2 Heisenberg model on the Kagome lattice, 
with both first- and second-neighbor interactions, is given 

by 

(ij) <fe» 
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FIG. S5: The entanglement entropy Sn(L y ) for the Kagome 



J1-J2 model in Eq. (|S3|) . with L y = 4 ~ 12 at L^ = 00. By 
fitting S n (Ly) = a n L y —"f n , we get 71 = 0.698(8) at J2 = 0.10, 
and 71 = 0.694(6) at J 2 = 0.15, while 72 « 0.42(1) at both 
J2 = 0.10 and 0.15. Here n is the entropy index. Inset: 



Here S* is the spin operator on site i, and (ij) ({{ij))) 
denotes the nearest neighbors (next nearest neighbors). 
In the numerical simulation, we set J\ — 1 as the unit of 
energy. We take the kagome lattice with periodic bound- 
ary conditions along a bond direction to define a cylinder, 
drawn vertically in the inset of Fig. IS51 and the unit of 
length equal to the nearest-neighbor distance. The re- 
sults for the entanglement entropy for J2 = 0.10 and 
0.15 are shown in Fig. IS5L for which both spin-spin and 
dimer-dimer correlation lengths are approximately one 
lattice spacing[9|]. We see that a linear fit using data for 
L y = 4 - 12 using Eq.© gives 7 = 0.698(8) at J 2 = 0.10 
and 7 = 0.694(6) at J2 = 0.15, both within one percent 
of In 2 = 0.693. In contrast, the extrapolated TEE 7„ 
based on the Renyi entropy (i.e., n > 2) clearly deviates 
from the expected value, although the cylinder's circum- 
ference is much larger than the correlation length, i.e., 
L y w 10£. For example, as shown in Fig |S5[ a linear fit 
using Eq.JTJ) gives us 72 « 0.42(1) at both J 2 = 0.10 and 
J 2 = 0.15, which differs from the expected value by an 
« 40% error. 



Kagome lattice with L x — 12 and L y 
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